Running-phase state in a Josephson washboard potential 
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We investigate the dynamics of the phase variable of an ideal underdamped Josephson junction in 
switching current experiments. These experiments have provided the first evidence for macroscopic 
quantum tunneling in large Josephson junctions and are currently used for state read-out of super- 
conducting qubits. We calculate the shape of the resulting macroscopic wavepacket and find that 
the propagation of the wavepacket long enough after a switching event leads to an average voltage 
increasing linearly with time. 



The dynamics of a large underdamped Josephson junction characterized by a capacitance C and Josephson energy 
Ej can be described by the motion of a particle in a washboard potential U{"f) = Ej{l — cos 7) -I- I^oj. The 
particle has C as the mass, the flux $o7 as the coordinate and the charge on the capacitor Q as the canonically 
conjugate momentum. Here 7 is the phase difference of the superconducting order parameter across the junction 
and ^0 = ^o/2tt = h/2e is the flux quanta divided by 27r. Much attention has been given since the discovery of the 
Josephson effect to the switching dynamics of the junction in the thermal activation regime and in the macroscopic 
tunneling (MQT) regime. Surprisingly, while the description of the state of the junction before a switching event 
and calculations of the corresponding probability has been a topical issue for many decades, what happens with the 
quantum state of the junction after tunneling did not receive that much attention. It is argued that the junction 
ends up in a running-phase state, with the voltage Q/C increasing until it becomes sufficiently large so that the 
transport could be done through quasiparticle excitations. However, a quantum mechanical description of this state 
is missing. Much of what we understand about the running-wave state, for instance the physics of the retrapping 
current, comes from assuming a quasi-classical dynamics. 

In this paper we give an explicit formula for the macroscopic wavefunction of an ideally underdamped junction after 
a MQT switching event. If the switching probability is exponential, which is the case for all the theoretical models 
and also confirmed experimentally, one expects (3| the following expression for the dynamics of the wavefunction 

l^'it)) = e-"/2e-'-»*|vI/o) + l^outit)). (1) 

In this equation, ^'o is the initial state, coresponding to a bound state inside one of the metastable wells, while ^'otjt(t) 
is the wavefunction of the particle corresponding to states in the continuum, outside the well (Fig. QJ. This expression 
gives indeed an exponentially decreasing probability for the particle to be inside the well, with lifetime F"^. In the 
following, we are interested in the structure of "^outit). 

To solve this problem, the standard approach is to start with a wavepacket localized initially in one of the metastable 
wells, and then expand it and evolve it in the eigenfunctions of the full Hamiltonian. This procedure works for simple 
potentials but even in these cases the solutions are complicated. Fortunately, unlike problems in scattering theory, 
in condensed matter the frequent situation is that we do not need an exact solution of the Scrodinger problem for 
tunneling, but rather we are interested in the most generic features of it. In most cases in solid state physics, tunneling 
is simply treated as a process that annihilates a particle on some mode of a solid and creates one on another mode. 
We will approach our problem in the same spirit j5| . A good approximation in MQT is that no other state within the 
well is involved with the exception of the state with energy loq in which the system is prepared, |5'o); therefore one 
can write a reduced Hamiltonian of the form 

i7-^c^o|*o)(*o|+ / ?ie|Ve)(V'.|+ / de[k{ujo,e)\^o){'>Pe\ + H^,LOo)\A){^o\], (2) 



where by {ipe} we denote the continuum of eigenvectors outside the barrier. We then write the wavefunction in the 
form 



with a(0) = 1, 5(0) — 0. Inserting this expression in the Schrodinger equation we get an integro-differential equation 
for a{t). The Laplace transform of this equation reads 



2 



where 

Xit) = ^ I de|fc(e,^o)|'e^(-«-)*. (5) 

In general, the tunneUng matrix element fc(wo, e) depend on the energies e and they are determined by the overlap of 
the left and right wavefunctions under the barrier ^J. We notice that since typically the lifetime of the metastable 
states is much larger than the oscillation period in the well (in other words the last term in the Hamiltonian is a 
perturbation), the states {|e)} which contribute effectively to tunneling are located in a relatively small energy interval 
compared to the plasma oscillation frequency, therefore the shape of these states under the barrier is approximately 
identical. We can then take the tunneling matrix element as being a complex constant; but since we will be interested 
exclusively in the outgoing component, the relative phase between the wavefunction inside the well and that outside 
will not play any role. We have confirmed this assumption also by expanding the initial wavefunction in terms of 
the WKB solution of the washboard potential calculated in Therefore we take k = hy^T/2TT real; we obtain 
£[x](s) = r/2= constant, which turns out to be the decay probability of the system. Indeed, the inverse Laplace 
transform of Eq. (^J gives precisely the classical exponential decay law 



The outgoing wavepacket becomes 



\^out{t)) = -^J- / de——— —li^t)- (7) 
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To conclude this derivation, we find that with the identification k — fi^Y jl-K the Hamiltonian Eq. Q becomes a 
model Hamiltonian for decay in the continuum which can be solved exactly, with solution given by Eqs. (^3) and Q. 
A similar type of model Hamiltonian has been obtained in 0,0. One can show, using the properties of the Lorentz 
distribution, that these wavefunctions are correctly normalized, as explained above. 
Let us now single out one component of the wave 'I'oMt(i), namely 
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We first notice that the normalization of the total function \'i>out) is such that (vPontl^oMt) = 1 — exp(— Fi), which 
reflects correctly the fact that the probability of finding the particle outside comes from an exponential decay law, 
while that of l^'^j) is such that (^'J^fl^'^J = 1. In the following we will see that l^'^j) plays indeed a special role. 
To move on, we notice that part of the expression for the outgoing phase contains a term which decays exponentially 




on a time scale F^^. These terms are associated with the fast components of the localized wavefunction which would 
escape first. Although a calculation that includes these terms is no doubt interesting, especially for the problem of 
non-exponential decay rates 4], in what follows we will regard them as transient oscillatory effects whose presence 
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will be difficult to assess experimentally anyway, and we will neglect their contribution. In the WKB approximation, 
far enough from the classical turning point, the eigenvalues {IV'e)} have the form (up to a normalization factor and 
constant phase factors due to matching to the region left of the classical turning point) of incoming and outgoing 
scattering states 
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(9) 



(10) 



The physical meaning of this voltage is that it corresponds to the (classical) energy accumulated on the capacitor when 
the phase difference across the junction is 7 and the initial energy of the system is e; indeed, CV^{'y)/2 = Tie — U{'^). 
We now use the fact that for values of 7 outside the well and far enough from the classical turning point the inequality 



\he — Uujq] <C Uujo — U{'j) holds. We can therefore take V^{-f) 
approximate the exponent as 



^0(7) ^(7) in the denominator of Eq. and 
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With these approximations, using Eqs. ((7|) and we can performe the integral over the angular frequencies e; as 
a result, the contribution of the in-going scattering states in zero, while the out-going scattering states build up a 
wavepacket of the form 
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Here A/" is a normalization factor which can be obtained through (i(i>o7)|^^t(7, i)P = 1 with the mention that 
we make a negligible error by extending the integral to —00, i.e. before the well region (where the actual values are 
exponentially small). The voltage ^^(7) is defined as ^(7) ~ Vb(7) — 'hLOQ/CVQ{'^) « Vb(7). It is interesting to see 
also what happens with the rest of the components of \^out)- Although they do contribute to the normalization as 
discussed before, they are decaying both in time and away from the barrier as exp[— VQ'^{ip)^Qdip\ which, as we 

will see below, would give far from the barier a factor of exp[— r/2(t -I- ^/2rj/ujp)]. It is clear that these terms can 
be neglected starting roughly from a time F"^. The wavefunction Eq. 1)12(1 contains all the information about the 
dynamical evolution of the state of the circuit containing the Josephson junction, and it is the main result of this 
paper. 

In a typical experiment, the voltage across the junction is monitored by a voltmeter at room temperature. A 
fundamental issue is to find a microscopic mechanism for the junction-voltmeter interaction and a suitable theory of 
quantum measurement that would model the collapse of the wavefunction; this is however beyond the scope of this 
paper. Still, Eqs. (|12|l and H13|l give a quite clear qualitative picture of what happens: the particle rolls down the 
washboard potential with a quasi-classical speed given by energy conservation C1q^(7)/2 = fiujQ — C/(7). Quantum 
mechanics enters in the picture through the tunneling rate; we expect the results of the measurements to have a 
spread dermined by F. One can assume that the measurement projects the outgoing state onto eigenvalues of the 
voltage operator; therefore the probability of recording the value V at the moment t will be given by the standard 
quantum mechanics recipe 



P{V,t) 
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As an example, had the outside the well potential U been zero, we would have gotten for the charge CV ^ by performing 
the integration in Eq. \i'6l , a standard Cauchy-Breit-Wigner distribution centered around CVq and full width at half 
maximum Th/Va 



P(CV) 



1 r?i 

n2Vo 



{CV-CVof 



Th 

2Vo 



(14) 



Considering again the case of a junction with a washboard potential U{'y), we notice that a good approximation 
is U{'-f) « Ej'j. This comes from the fact that switching is typically observed at values of the bias current close to 
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the critical current of the junction, as well as from the observation that for times larger than T~ the wavepacket is 
concentrated at large values of 7 ^ 1, in which case the cos 7 term in the potential is negligible. In other words, the 
particle gets soon so fast that the " speed bumps" created by the Josephson effect are not slowing it down significantly. 
This can be checked a posteriori. A first observation is that the relevant quantity for the dynamics of the center of 
the wavepacket is the argument of the O function; the condition that this argument vanishes sets the maximum value 
of l^'outp and gives a phase LUpt/2 ^ 1 for i larger than F^^. A legitimate concern is whether the wavefunction does 
not spread faster than it moves downwards. This is not the case, as we will see below: the spread of the wavefunction 
increases linearly with time, while the average coordinate (phase) is advancing as t^. With these observations, the 
normalization constant can be calculated, and the outgoing wavepacket becomes 



e (tLUp - . (15) 

A plot of the wavefunction is given in Fig. |21 
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FIG. 2: Contourplot (7, t) of the modulus of the running-phase wavefunction Eq. for cjp = 60r — 30 GHz. 

The average phase (flux) corresponding to this wavepacket can be obtained 

(7)(i)=a;^(tV2-r-4 + F-2), (16) 

and we notice that the dominant term is quadratic in t. The spread of the flux variable is given by (we keep only the 
dominant term here) 

V{l'){t)~{7m^u;lr-H. (17) 
To get the average voltage we can use Ehrcnfcst theorem; we obtain 

{V{t)) = ^,ujl{t-r-^). (18) 

The dominant term for the voltage is linear in time and satisfies the classical energy conservation C{V)'^{t)/2 = 
n^o-C/((7)(i))«/o(7)(0- 

Let us now analyze what happens in typical switching current experiments, as they are done now in the context 
of superconducting qubits 0: the bias current of the junction is increased fast to a value that allows tunneling, it is 
kept there for a time < r < F~^, then it is lowered to a value that suppresses tunneling. This value has to be large 
enough so that the experimentalist can get a reliable reading of voltage on the quasiparticle branch if the junction 
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has switched; in practice, it can still satisfy / « Jq. Although the change of the bias current has a major effect with 
respect to tunneling through the barrier, where the the tunneling rate decreases exponentially with the height of the 
barrier, from the point of view of the structure of the running-phase state it amounts only to a modification of the 
parameter F. Finally, the current is put to zero and, after waiting long enough for retrapping to occur, the whole 
cycle can be repeated. In our model, the essential physics is that after the time t, the tunneling matrix element t is 
zero, therefore the system evolves only under the action of Hq. The wavefunction is "cut" into two separate pieces, 
one which is (almost) the bound state ^'o inside the well, the other being the wavepacket in the continuum which 
evolves as 

/ r r p(-T /2+ibJo)T -ibjj(t-T) _ p-iet 

l^outit)} = - / de — — ^ (19) 

\ 27r J F/2 + i[ujo - ej 

with normalization (^oMt(Ol^o«t(0) — 1 ~ e^^'^. Now, for t — t > F^^ we can see that the outgoing function consists 
of two consecutive (separated by the time t) and dephased (with loot) outgoing wavepackets with the structure of 
l^ont(^)) which propagate at the same speed across the phase coordinate 7. The second wavepacket, which has a 
probability amplitude smaller by a factor of e~^/^'^, results from the waves localized near the barrier during the time 
T when tunneling was in progress. After integration over energy, we get 

\^,^,{t)) = e-(r/2+-«)n*rn.(i - r)) \^7ut{t)). (20) 

where \^7ut{^)) given by Eqs. ifT^ and ifT^ . To check that the normalization {'^ out{t)\^ out{t)) 1 — e^^'^ 
remains valid, we notice that in the region of overlap of the two wavepackets, which coincides with the domain 
where |\['^f(t — t)) is finite t — t > ^QV7^{ip)d(p, there exists a very simple relation between them: |^'^t(i)) = 
exp[— jwoT — FT/2]|'J'^((t — r)). Using this property and the previous expressions Eqs. lfTC|) and ((TH)l . we can calculate 
the average phase and voltage on the state Eq. H20() 

(7)(t) = - e-^^]e - e-^^u^lrt - ^e^^^^^, (21) 

and 

{V) {t) = $oc^pt (1 - e-^^) - ^^ule'^W. (22) 

In Fig. 131 we present a plot of the average voltage as a function of the time r. We see that for values of r of the 
same order or larger than the lifetime F~^ the average voltage at t flattens, reflecting the fact that the junction has 
switched, as in the case of Eq. H18() . 
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FIG. 3: The average voltage as a function of the time r for — 30 GHz, F = 1 ns and t = 10 ns. 



For designing an experiment to test these predictions, several remarks should be made. In the case of real junctions, 
the Josephson energy and the plasma frequency can be reduced by using a SQUID configuration and by adding 
capacitors in parallel with the junctions. This makes the time evolution of the switching state slower and therefore 
easier to detect. An important limitation on time comes from the fact that as soon as the voltage reaches the 
quasiparticle branch (at twice the value of the gap) our analysis is not valid. The other limitation is technological: 
even with a good dilution refrigerator, thermalizing the junction is very difficult at low temperatures. With a good 
high-power refrigerator with base temperature of about 5 niK, we assume an optimistic value of 10 mK for the effective 
temperature of the electrons. This temperature corresponds to a crossover angular frequency of 8.66 GHz between 
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the MQT and the thermal activation transition. A plasma frequency of LOp = 30 GHz (zero bias current) will thus 
keep us safely in the MQT regime when the current is raised up to about half a percent close to the critical current, 
according to the formula that gives the plasma oscillation frequency at a finite bias current P, Q ■ For Nb, with gap of 
1.4 meV, this corresponds to a time of approximately 10 ns, as given by Eq. (|18|) . A voltage increase on this timescale 
can be detected with standard experimental techniques Suppose now that we choose to work at currents about 5% 
less than the critical current. We still have to satisfy the condition t > F^^; an inspection of the formula that gives 
the tunneling rate for underdamped junctions (see e.g. 01) shows that switching rates of about 500 MHz and more 
(with the restriction F <^ cUp) can be achieved for Ej/tUp of the order of 30, values which can be obtained easily with 
large junctions. 
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